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Abstract: We study the dynamics of light-like tachyon condensation in a linear dila- 
ton background using level-truncated open string field theory. The equations of motion 
are found to be delay differential equations. This observation allows us to employ well- 
established mathematical methods that we briefly review. At level zero, the equation of 
motion is of the so-called retarded type and a solution can be found very efficiently, even 
in the far light-cone future. At levels higher than zero however, the equations are not of 
the retarded type. We show that this implies the existence of exponentially growing modes 
in the non-perturbative vacuum, possibly rendering light-like rolling unstable. However, a 
brute force calculation using exponential series suggests that for the particular initial con- 
dition of the tachyon sitting in the false vacuum in the infinite light-cone past, the rolling 
is unaffected by the unstable modes and still converges to the non-perturbative vacuum, in 
agreement with the solution of Hellerman and Schnabl. Finally, we show that the growing 
modes introduce non-locality mixing present with future, and we are led to conjecture that 
in the infinite level limit, the non-locality in a light-like linear dilaton background is a 
discrete version of the smearing non-locality found in covariant open string field theory in 
flat space. 
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1. Introduction 

While the free action of open string field theory Q is local in the sense that it involves not 
more than two derivatives of the string field, it is well known that the interaction term of 
Witten's string field theory [Q] contains infinitely many derivatives. Theories with more 
than two, but finitely many, derivatives []|, || suffer either from an unbounded Hamiltonian 
or from ghosts 1 , but these instabilities do not necessarily survive in the limit of infinitely 
many derivatives. Perhaps the simplest way to see this is that the propagator in a theory 
with finitely many derivatives is the inverse of a polynomial and therefore has poles, some 
of them ghosts. In the limit of infinite number of derivatives, however, the propagator 
becomes the inverse of a function that might have only one zero, corresponding to a reg- 
ular excitation. It might even have no zero at all, like in p-adic string theory where the 
propagator is an exponential (furthermore, this exponential propagator renders all loop 
diagrams finite ||). 

In cubic open string field theory, the form of the nonlocality is universal, the higher 
derivatives of any field <p{x) always appearing in the interaction as 

0(x) = ET n 0(x), (1.1) 

1 See, however, JBJ for an example of quantization of a theory with four derivatives, with no ghost and 
bounded (but non-hermitian) Hamiltonian. 
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where K = ^p, and we use the signature r}^ v = diag(— 1, 1, . . . , 1). In this case one can 
see the nonlocality explicitly because <j>{x) is a smearing of (j){x) as can be seen from the 
convolution formula J5J 

e^(j)(x) = —= / e ~^ x ~ v) <f>(y)dy, (3 > 0. (1.2) 



For a homogeneous time-dependent problem, one would take <j)(t) as the fundamental field 
and write <f>(t) = e lo ^ 9 ?4>(t), which can then be written as a convolution as in Eq. §L%). 
A consequence is that the equation of motion of a homogeneous time-dependent string 
field involves the string field not only at time t but at all times, both in the past and in 
the future of t. One can treat this kind of equation either as a convolution equation using 
Eq. (|1.2| ), or as a differential equation of infinite order. It must be understood, however, that 
such a differential equation cannot be seen as a limit of finite-order differential equation; 
in particular the initial value problem becomes different when we have infinitely many 
derivatives (see || for a rigorous discussion, and || which contains some similar results). 

A time-dependent equation of motion with infinitely many derivatives which is of 
particular physical interest, is the equation describing the decay of an unstable D-brane. A 
well-known problem is that, on the one hand, a boundary conformal field theory (BCFT) 
analysis shows that one should expect a monotonic decay of the tachyon down its potential 



|1C] and that the energy of the D-brane is converted into very massive closed strings at 
rest, behaving like dust |ll[] (tachyon matter). On the other hand, numerical solutions of 
string field theory |l^, [D| show a completely different behavior. Namely, the tachyon does 
reach the non-perturbative vacuum, but it then continues further and starts oscillating 
around it with diverging amplitude. Although the tachyon can climb arbitrarily high up 
the potential, the energy conservation is not violated because the kinetic energy can be 
negative. In fact, it may seem that energy conservation obviously discards a monotonically 
rolling tachyon that would stop at the local minimum of its potential. However, this is not 
totally trivial because one could imagine that the energy of the D-brane is somehow stored 
in the very high-order derivatives of the tachyon, still allowing for a monotonic rolling. 



But this is actually ruled out |14j, |l5j because one can write the expression of the energy 
in an integral form which makes it clear that, in fact, a monotonically rolling tachyon 
cannot conserve energy. It is now believed that these ever-growing oscillations are not 
catastrophic after all. For one thing the string field is not a gauge-invariant observable; 



but more concretely, it was shown in [16], that a field redefinition mapping the cubic SFT 



action to the boundary SFT action, would also map the oscillating solution to a well- 



behaved solution. More recently, it was shown [17] that the closed string boundary state 
obtained from the rolling tachyon solution, coincides with the BCFT boundary state. 

It is interesting to investigate how this wild rolling changes if we somehow couple 
the closed strings sector to the open string SFT action. The most consistent way to do 



this would be to consider open-closed string field theory [18]; but solving the equations 



of motion of the purely closed sector |[{j] involves a much higher level of difficulty [20 



21, 22, 23, 24, 25, [26|, ^7], 28]. A somewhat more manageable approach would be to 



consider a fixed closed string background, but the SFT action becomes in general non- 



- 2 - 



polynomial in a generic closed background R29|| , hence also hard to solve. What can be 
done, however, is to minimally couple gravity to the open SFT action. It has been shown, 
for instance, that minimally coupling an open superstring tachyon to a FRW metric tames 
the wild oscillations of the tachyon; and convergent rolling tachyon solutions were found 
numerically JO]]. In |3C], which was the motivation for the present work, Hellerman and 



Schnabl considered open SFT in a linear dilaton background. They chose a light-like dilaton 
gradient and a string field depending only on the light-cone time x + . This is physically 
motivated because a bubble of true vacuum is expected to expand at the speed of light |H| . 
If the radius of the bubble is large enough, we can focus on one small patch and approximate 
it by a plane, and we choose the light-light coordinate x + (which we call light-cone time) 
to be orthogonal to this plane. Moreover, with this ansatz important simplifications occur. 
In particular, using the fact that e x+ is an exactly marginal operator, Hellerman and 



Schnabl were able to use the results of [32, |33| in order to prove that the rolling tachyon 
asymptotes to the tachyon vacuum at large light-cone time. On a more explicit footing, 
Hellerman and Schnabl also considered the SFT action truncated at level zero in Siegel 
gauge (i.e. keeping only the tachyon). Here the light-cone simplification manifests itself 
by changing the nature of the non-locality. The non-locality of Eq. (1.1), which by virtue 



of Eq. (F2), involves the tachyon field at all times, becomes simply the tachyon at some 



retarded light-cone time <p(x + — 7). The equation of motion for the tachyon is then 

(j)'(x + ) -cp{x + ) = -K 3 0(x + - 7 ) 2 . (1.3) 

A numerical solution to this equation was worked out by Hellerman and Schnabl. Their 
method, however, didn't allow them to go very far in light-cone time, but enough to 
see convincingly that the tachyon reaches the vacuum after oscillating around it with a 
decreasing amplitude. In f35|| , Barnaby et al. considered the initial value problem and the 
stability of light-like rolling in p-adic string theory and in SFT at level zero. In particular, 
they were able to numerically solve Eq. (|1.3|) for a much larger light-cone time interval. 



Their numerical method is based on the diffusion equation |3^, 37, |38| . Although this 
method can in principle be generalized to higher levels, it is hard to do so in practice. 

The motivation for our present paper, was to investigate further the light-cone rolling 
in Siegel gauge by considering higher-level fields. We will consider levels (2,4), (2,6), 
and (4,8), where the notation (L, M) means that we are keeping fields up to level L and 
interactions up to total level M. Our results are three-fold. 

Firstly, we realized that equations of the type (|1.3j ) are known in the mathematics lit- 
erature as delay differential equations (abbreviated DDEs). The most widely used method 
for numerically solving DDEs is the method of steps. Using this method, we show that 
solving Eq. (|ll^) numerically becomes surprisingly easy. Moreover, the generalisation to 
higher levels is straightforward. 

Secondly, we show that when we include higher- level fields, the nice picture of the string 
field gently oscillating around the vacuum with decreasing amplitude, takes a serious hit. 
Indeed, we show that already at level two, the tensor fields which must be included in our 
analysis, bring derivatives into the equations of motion in such a way that these become a 
so called system of higher order neutral DDEs with one positive delay Such DDEs cannot 
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in general be solved with the method of steps. We have to do some simplifications before 
obtaining numerical solutions. What we can do, however, is to look at the equations of 
motion close to the vacuum. We will find that the latter effectively contain several delays. 
This would pose no further conceptual difficulty if all delays were positive, but we show 
that we obtain negative delays as well. In other words, the equations of motion effectively 
involve the fields at some past light-cone times, but also at some future light-cone times. 
We show that if we have both negative and positive delays, there exist growing oscillation 
modes around the vacuum. This suggests that the string field may not converge to the 
non-perturbative vacuum. This seems to be in contradiction with the analytic solution 
obtained in [30], and also with our numerical solution obtained by expanding the fields 
in exponential series (this provides in principle a very accurate solution but only up to a 
limited light-cone time). But the two pictures can be reconciled if we notice that the initial 
conditions for the analytic solution (which are the same as those for the exponential series 
solution) are very special. The diverging modes might not be excited for this particular 
solution, but our results imply that a small change in the initial conditions can render the 
rolling non-convergent. 

Thirdly, by studying the equations of motion near the vacuum at level four, we show 
that there are more delays at this levels, and that they are more spread, both towards 
the past and towards the future. We are led to conjecture that in the large level limit, 



one recovers a discrete version of the non- locality (1.1). Subsequently setting the dilaton 



gradient to zero, the delays become less and less spaced, and we will recover (1.1). We 
conclude that in particular, the nice simplification of the non-locality that happens at level 
zero, is only accidental. 

This paper is structured as follows: In the next section, we calculate the action and 
derive the equations of motion at level two. We give a short review of DDEs in Appendix 
^ and the complete results are given in Appendices |B] and O. We solve the level-zero 
equation with the method of steps, and we attempt to do the same at level two. We show 
what problems we face and what can be done to get some information on the rolling at 
this level. In Section |3], we study the linearized equations of motion near the vacuum. We 
show that negative delays appear at level two and four and conclude the general form of 
non-locality in the linear dilaton background. In order to do so, we need to calculate the 
determinants of polynomial matrices. A naive approach fails if the matrices are too large, 
so we explain a little-known method for calculating such determinants in Appendix In 
Section ||, we discuss further the consequences of our results. And at last, we offer a review 
of linear dilaton CFT in Appendix ||. 



2. Light-like dynamics of the vacuum transition in Siegel gauge 

This section is divided into several paragraphs. At first we show how to derive the equations 
of motion of open string field theory in a linear dilaton background with level truncation. 
We present more details of the calculation, on the one hand to introduce the notation, 
and on the other hand because they are omitted too often. After that we briefly describe 
how we solved the resulting delay differential equations on the computer. To that end we 
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explain how astonishingly natural it is in our setup to choose the infinitely many initial 
conditions required for producing a unique solution. We test our machinery in the simplest 
possible case of level zero, and observe excellent agreement with the literature |3^, |35| . 
We then go beyond level zero and explain our results at level two and four, which can 
be summarized as follows: The individual modes of the string field are initially in the 
perturbative vacuum, then, driven by the tachyon, they grow steeply. Finally they oscillate 
around their respective vacuum expectation values with decaying amplitudes. 

Derivation If we write the string field in terms of vertex operators as = \I/(0)|0), the 
action of Witten's open string field theory reads 

S= ~7 Q(^W + ^(/l°*( )/2 O *( )/3 o *(0))) , (2-1) 

where the functions fi are the conformal transformations mapping each string (semi-disk) 
to the common interaction upper-half plane. 

We have derived the action and equations of motion with two independent methods. 
First we calculated the action by hand, calculating explicitly the conformal transformations 
fi o \&(0), and then the CFT correlators 2 3 : 



v 



Yi :e *-*c«) : jj a^x* (4) ) = (2nf ^ n 

\i=l j=l I \ i / i„j=l,i<j 



1 2 cx /c ^ ■ kj 
i 3 1 



•^n^(4)+^(4)] / >- (2.2) 

The new objects v, q serve as a tool to quickly work out the combinatorics of the con- 
tractions - just expand the product into a polynomial in v, q and observe the following 
rules: 



1. Replace {£) = —iol 5^=1 S^F- 

2. Contract products of two qs using —a'r) fJ,l '/2 (z — z') 

3. Remove all terms with an odd number of qs. 



4. Note that the general expression diverges if z = z' . For correlators of normal ordered 
products (e.g. *QX fJ, e lk ' x l x *...*) these terms precisely cancel, providing a non- 
singular result. In this case we can further simplify with the additional rule 

Remove all terms in v with z = zt and all q-products at the same point, z = z' . 



The dilaton background enters explicitly in two ways: 

2 For a detailed derivation of the correlator in the linear dilaton background, cf. |39] 
3 Note that we use the complex derivative d z instead of the real derivative. 
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1. the delta function is updated to include the breaking of the translation invariance by 
the linear dilaton background 




where the delta function of a complex argument is formally defined by the following 
integral representation 

S D (^ki + iv) ee ^p/d D . v ". (2.3) 

2. through the modified conformal transformation law 

X» (z, ^/ol" (z, z) = X» (/ (z) , f (z)) + log \f (z)\ 2 

needed when mapping the string field vertices to the interaction worldsheet. 

In order to make sure that our results are correct, we redid the same calculation with 
the method of conservation laws. Luckily, the conservation laws for an anomalous vector 
(like dX^ in a linear dilaton background) were already worked out by Rastelli and Zwiebach 
in p5| . This method has the advantage of being easy to implement on a computer. In our 
case we wrote a mathematica program. To our satisfaction both methods agreed entirely. 
Since the calculations become rather cumbersome at higher levels to do by hand, we relied 
on our code for the equations of motion at levels (2,6) and (4,8). 

Explaining our notation we will quickly see that at level two, the string field written 
out in its mode expansion can be reduced to just eight spacetime fields. We follow the 
conventions by pl| (except that we call (3 their f3\) and write the string truncated to level 
two fields. Working in the Siegel gauge, we can eliminate those terms containing a cq -ghost 
mode, and due to the twist symmetry of the action, we can consistently set all terms at 
odd levels to zero. This leaves us with the following expression for the string field at level 
two: 

|*) = j</>+ -^BfM&-2 + -^ B ^ a> ^i a ~i +/36_ic_i|ci|0). (2.4) 
Since we work in D = 26 spacetime dimensions, expr. (2.4) contains 379 spacetime 



fields. But in the case of light-like tachyon rolling, we can drastically reduce the number of 
fields needed in our calculation. Working in the light-cone frame we can split the dimensions 
into light-like and ordinary components 

fi= (0,1,2,... D-l)^ (+,-,2,3,... D-l) = (+,-,*). 

We assume the linear dilaton gradient light-like, V 2 = 0. By rotational symmetry we 
can choose a coordinate system where V = (V + , 0, . . . , 0). Furthermore we consider space- 
time fields 4>, (3, B^, B^ v in expr. ([O]) that depend only on the first lightcone coordinate 
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x + . As detailed in appendix |B|, we can then focus on the following eight (of 379) fields to 
compute the action: 

{<f>,B + ,B-,B++,B + -,B—,F,p}. (2.5) 

Note that is the tachyon field and F is the scalar field associated with the contribution 
of BV, i,j = 2. . . 25 to the trace of B^ by 

Tr B fJ,v = —2B + ~ + F. (2.6) 

The resulting action is presented in appendix [B[ both in Lorentz covariant form and ex- 
plicitly using fl2.5|) . To check its correctness, one can take the limit of vanishing dilaton 



gradient, V — )■ and compare to the action found in [41|. Both expressions agree as desired. 

With the action computed, we can proceed to deriving the equations of motion in 
notationally compact manner. For a Lagrangian £(</>, d(f>, d 2 (/)... ) containing arbitrary 
orders of field derivatives d n (f), the Euler-Lagrange equation is 

dC dC dC 

= — " fyi a ra IT + 5 mA 2 



In this notation the derivatives are not symmetrized, their order matters: 

d[d Vl d U2 ...d Uk 4>] w w 

This is just a matter of more convenient bookkeeping as we sum over all combinations of 
indices. In a compact notation we can define the differential operator T>^ which returns 
the equation of motion for the field (f> when applied to the Lagrangian C depending on (j> 
and possibly other fields to any order in the field derivatives: 

vt - g ( - 1)t9A - 9 ^[^...^] 

V^C = 0. 



Let's see how to apply this in a concrete example, take a generic interaction term from C, 
e.g. A(x + )(d^ 1 . . . <9 W B(x + ))(d^ l . . . d fll (f)(x + ))e v+x and compute its contribution to the 
equation of motion. When deriving the equation of motion for the tachyon field (j> we apply 
T)® . We use □ = —2d + d- + d{d l = —2d + d- because £ is independent of the Xj-coordinate. 
Fields with a tilde are defined by 

oo 1 

4>(x + ) =K a ' U (f ) {x + ) = Y,- } {n'^gK) n {-2d + d_) n (t>. 

n=0 n ' 
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Now derive the equation of motion: 

00 1 

= {-l) k d Vl ...d Uk -{-2a'\ogK) n 5^...5^...5^5^5 2n+l)k 

k,n=0 
00 1 

= E(-!) 2n+ ^ Ha'iogK)"^ ...a tt (a+a_) B 

n=0 

= (-i) l ^...a w e- 3ta ' to «(w-. 

Applying 2?^ to the whole term, it becomes apparent that it acts essentially as a translation 
operator when cL — > V + : 

V* [A(x + ) (> . . . d^B{x + )) (d^ . . . e v+x ~ } 

= (-l)'a^ ...3 w e- 2a ' log W 9 + a - (V 1 ...d w 5(x+)) e y+x "} 

= (-1)' W . . . 9 w e- 2a, ^ + W+ (d w . . . d^B{x + )) e v+x ~ } 

= (-1)' M1 . . . W {A (x + - 2a'y + log (if)) (9^ . . . d w £ (x+ - 2a' V + log (if))) e l/+x " } 

where finally the tilde was eliminated because DA(x + ) = 0. The translation operator 
simply shifts the argument x + by — 2a' V + log (K). This is a generic feature in every 
interaction term, we therefore define the symbol 

y + = x + - 2a'V + log (K) 

for the shifted point to abbreviate the notation. The reason for introducing T>^ is that it 
simplifies the calculation significantly. As an example, consider the chain rule 

V* [A (x + ) 2 (x+) e v+x ~) = 2A <p (y+) e v+x ~ . 

The set of eight equations of motion contains derivatives of up to fourth order, each equation 
is quite lengthy with one notable exception: the equation of motion for B is short and 
can be solved easily: 

V*-L total = + B++'(x+) + * B++ ^y + ) 1 0. 

The equation is linear in B ++ , hence we set B ++ = to obtain a solution. As a matter of 
fact this still admits non-trivial solutions for the other seven fields. This observation was 
also made in p^], in the related context of OSFT using the lightcone basis for the modes 
of the string field. This is a peculiarity of level (2,4). Indeed, already at level (2,6) one 
cannot consistently set B ++ to zero anymore. 

Setting B ++ to zero at level (2,4) reduces the length of the equations of motion by 
about one third and the total differential order from 21 to 17. The resulting set of equations 
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of motion is presented in appendix |C|. The seven equations contain a total of 144 terms. 
For reference we list the seven remaining fields 

{0> B + ~, B , F, /3, B + , B~, } . (2.7) 

Let tyj, j = 1 ... 7, denote one of the above fields, e.g. 9i (x + ) = 4> ( x+ ) is the tachyon 
field. Then the generic form of an equation of motion is 4 : 



= (x+) + 89 j (x+) + E E a nm^i (y + ) d m 9 k {y + ) . (2.8) 

i,k=l n,m=0 

The first two terms come from the kinetic term in the action. They have the same sign 
except for the tachyon and are evaluated at position x + . All derivatives are understood 
with respect to x + . All terms in the sum arise from the cubic interaction. They are 
evaluated at y + and contain derivatives of up to the third order. In fact many of the 
real-valued coefficients a l ^ m are zero. 

Solving the Equations of Motion Given the general structure of the equations of 
motion fl2.8| ), let us now focus on solving them at level (2,4). There are derivatives with 
respect to x + at both points: x + (max: 1st order) and y + = x + — 7 (max: 3rd order). 
The total differential order of the system of equations is 17. Because of the non- locality 
(fields at x + and at y + ) these are not ordinary differential equations (ODEs). It is the key 
observation that, mathematically speaking, we have a system of coupled delay differential 
equations of the neutral type (NDDE) with one constant delay 7. The theory of delay 
differential equations (DDE) has been developed to great extent in the 20th century as this 
type of differential equation arises in a large array of disciplines: populations dynamics, 
machine control theory, neutron diffusion, spreading of diseases, retarded propagation in 
classical electrodynamics etc. We give a short introduction to DDEs, highlighting the 
differences to ODEs, in appendix [A]. Further useful references are [43, 44, [45| ]. 



In order to solve DDEs we use the method of steps, explained further in Appendix [A|. 
The basic idea is to reduce the problem of computing the solution over a full interval to 
subintervals where the DDE reduces to an ODE which is solved with standard methods. 
In order to obtain a unique solution, it is however not sufficient to give an initial condition 
at one point as in the ODE case: initial conditions over a finite interval (of the length of 
the delay 7) have to be specified. These conditions are called the initial data. They are 
essentially uniquely determined when we require that the tachyon condensation starts in 
the perturbative vacuum. 

There are other methods. For instance Barnaby et al. transformed the level- 

zero DDE into an equivalent diffusion-like local partial differential equation problem with 
suitable boundary conditions, solved that with standard codes and finally transformed back 
to obtain the DDE solution. For the system of equations at level (2,4), their method seems 
cumbersome to apply. In contrast, the method of steps is powerful enough to easily extend 



4 To simplify the notation, we set V + = a' = 1. 
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to the case of several unknowns. For more details on diffusion methods, we refer the reader 
to 

We have seven coupled equations for seven fields with derivatives up to the third 
order. The numerical stability and convergence properties of the method of steps have 



been studied carefully in the past decades, cf. [45] for an extensive review of numerical 
methods. The problem can now be considered a standard one, much like solving a system 
of ODEs numerically. This implies that the user has to carefully check the consistency 
of the numerical solution, preferably by independent methods. Much of the effort in this 
section is aimed in that direction. 

On the computer we need dimensionless numbers. As of now we work in units where 



a' = 1. For convenience we choose V + = 1. 

Choosing initial data Let us now see how to supply initial data for each field to obtain 
a unique solution to the equations of motion. A priori we could choose any initial data, 
but physical reasons nearly completely fix them. Consider as an example the tachyon field 
(j). In principle we would like to obtain the solution <p{x + ) Vx + 6 M. from the equations of 
motion. On the computer we can only obtain a solution on some finite interval [a^ in , £+ax] • 
We then need to fix the tachyon on [a^m ~~ 7' x miJ' w hh the delay 7 = 2 log if. To find 
constraints recall that we are interested in the tachyon condensation solution: the solution 
should initially be in the perturbative vacuum (string field |\&) =0) and in the end arrive 
at the non-perturbative vacuum. Hence at sufficiently small x^ in the absolute value of the 



tachyon and all other six fields $/j(x + ) from (2/7) should be small 



As we take a^in ~~ ^ ~°° we can neglect all interaction terms ~ SS?i&j in the equations of 
motion (|2.8| ) and receive dramatically simplified equations. In other words, we linearize 
the equations of motion around zero, = + , and neglect all terms quadratic in the 
small quantity S^fj. The resulting equations are 

= 5Vj (x + ) + (-l)^* 1 dSVj (x + ) . (2.9) 

These are first order ordinary differential equations which can be solved analytically, we 
just need to fix the initial condition. Starting in the unstable vacuum means all fields 
vanish at negative infinity, 

*j (-00) = 0, Vj. 
The linearized tachyon equation of motion has the solution 

5(f> (x + ) = a 6 e x+ (2.10) 
= 5(j) (—00) =>• arbitrary. 

The solution grows exponentially, reflecting the instability of the perturbative vacuum. The 
initial condition does not fix the constant as- This is the only free parameter in choosing 
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the initial data as we will see shortly. In short, there are just three cases, 



a<j> 



positive 


negative 

that give qualitatively different solutions to the full equations of motion. From equation 
( 2.10 ) gu can be chosen to have any real value. If we set it to zero this corresponds to the 



static solution in the perturbative vacuum that we are not interested in. 

For 7^ we obtain interesting solutions. By shifting the origin in the x + direction, 
we can fix the absolute value 

M = 1 (2.11) 

without loss of generality. It turns out that for < the solutions diverge. This can 
be intuitively explained as "rolling down the wrong side of the hill" : the effective tachyon 
potential is presented schematically in Fig. p]. The linearized equations of motion for the 



V{<t>) 




Figure 1: Effective tachyon potential. Depending on the sign of the coefficient the tachyon either rolls 
off to infinity or to the non-perturbative vacuum. 



other fields have the solution 

8Vj (x+) = aj e~ x+ 

= (-OO) => GLj = 0. 

In conclusion the initial data are as follows: the tachyon rises exponentially with a prefactor 
of choice, all other fields vanish. This confirms that the tachyon drives the condensation 
process. The same in a formula is 

5Vj (s+) = 5 jAa<p e x+ , x+ n - 7 < x + < x+ n . 

We require that the solutions of the equations of motion be analytic, thus we can express 
the analytic initial data as 
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Method 


AccuracyGoal 


WorkingPrecision 


MaxSteps 


Adams 


15 


20 


50000 



Table 1: Standard options used with mathematica 7's NDSolve to numerically solve the equations of 
motion. 



With our choice for the tachyon S^i = a^e x , we have fixed every coefficient, c\ 
Conversely for the other fields we find = 0. 

We conclude: the solution of the linearized equation requires only one initial condition 
(a<£ for the tachyon) , it then fixes the countably many initial conditions required for finding 
a unique solution of the full non-linear DDE. Note however that in general a Lagrangian 
with derivatives of all orders does not require supplying countably initial conditions for a 
unique solution, cf. Q for a review of the initial value problem for linear equations. 

Programming details Several codes for NDDEs implementing the method of steps are 
available. We choose to use mathematica in the version 7 as it easily allows to further 
manipulate the equations symbolically besides the capability of solving systems of NDDEs 
with the single command NDSolve. 5 mathematica also allows to do the numerics with 
arbitrary precision, which we made use of, as the built in machine precision was not quite 
satisfactory. We want to warn the reader that in the case of NDDEs with higher order 
derivatives at delayed positions mathematica quickly returns results without any warning 
or error message. However upon plugging the supposed solutions into the equations of 
motion we realized that they do not satisfy the equations. As always when using numerical 
results, checking is crucial. In those cases where mathematica yields correct results we 
used the parameters and options listed in Table [j]. 

Warm up at level zero As a basic consistency check for our numerical method we 
run the simplest example: the level zero truncation to the tachyon only. We can compare 



the results to Barnaby et al. |35| (diffusion problem) and Hellerman and Schnabl [30] 
(exponential series solution) that each solved the same problem with a different method. 
The equation of motion |}(]] for the tachyon at level zero is 

= (f) 1 - <f) (x + ) + K 3 (j) 2 (y+) . (2.12) 



We used the initial data 0(x + ) = 1 • e x+ , Eq. ( 2.10| ), on an initial interval [—25 



2 log X, —25]. In practice this is close enough to the perturbative vacuum, as (p{x + ) = 
e -25 ^ ^ Moreover, the number of intervals between x + = —25 and x + = is 2 v^,k ~ 
47.8. This tells us (see Appendix [A|) that the numerical solution will be differentiable at 
least 47 times for x + > 0; we can therefore expect that it will be a very good approximation 
to the analytic solution. The solution is depicted in Fig. [2[ It is in excellent agreement 
with solutions from p0| , p5| . The asymptotic behavior for large x + is known, it is of the 
form exp x cos. The frequency and decay rate can be determined by a simple fit, the 
values agree perfectly with those noted by Hellerman/Schnabl from linearization at large 



3 This feature was not available in previous versions. 
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Figure 2: Level zero tachyon condensation calculated with the method of steps. From the old vacuum 
<j> = we jump to the new vacuum with exponentially dampened oscillations. The vev is indicated by the 
dashed line. 



x + . Note that these authors could compute the solution of ( 2.1 2| ) only as far as x^ ax = 7 
because of computing time limitations: the computational complexity grows exponentially 
with x + for their method. Even though we plot the solution only up to = 25 we 

compute it up to cc^ax = 100 and beyond without any difficulty The numerical solution 
(~0.1s) actually takes less computing time than plotting the result (~ls) on a modern 
computer. This confirms that the method of steps works both fast and accurately. 

We want to emphasize that the level zero equation of motion is in many facets simpler 
than the level two equations of motion. Obviously it is only one equation compared to seven 
coupled equations, but the major difference is of another kind: the level two equations are 
neutral DDEs, while ( 2.12j ) is of the retarded type, there are no derivatives of (/> at position 

y + . 

Level two solutions The full set of seven equations of motion at level (2,4) (appendix |(|) 
contains derivatives up to the third order in the fields at the delayed position y + = x + — 7, 
and up to first order at x + , the total differential order is 17, while the total number of terms 
is a staggering 144. As explained above mathematica 7 cannot handle higher derivatives 
at delayed positions, and as of the time of writing we know of no other numerical method 
for solving a system of higher order neutral DDEs. Hence we look for ways to simplify the 
problem that allow us to follow the vacuum transition: 

1. Set all higher field derivatives to zero. 

2. Consider only the scalar fields, <j),/3,F, and set the other vector/tensor component 
fields, , B~~ , B + , B to zero. Then the system of equations contains at most 
first order derivatives. 



-13- 




1st Deriv. 
4> Scalar 
4> Exp. series 

(<!>) 



Figure 3: The different approximations for the tachyon condensation at level (2,4): From the perturbative 
vacuum <f> = the tachyon jumps to the non-perturbative vacuum. The vev of the full system is indicated 
by the dashed line. When ignoring higher derivatives the vev is unchanged. In contrast when considering 
the reduced system composed of only the scalar fields <j>, F, fi the tachyon vev is slightly altered; this can 



be understood from the remark starting before Eq. (2.14) 



3. Rewrite the fields as exponential series, 

oo 

*;(* + ) = X>> em:+ > ( 2 - 13 ) 

71=1 

and solve for the first few hundred coefficients a,j^ n recursively. Evidently for large x + 
this procedure requires knowing many of the Oj jn , in fact the number of coefficients 
needed for an accurate solution grows exponentially with x + . Thus we compute 
the solution with this approach in reasonable time only on a relatively small range 
[^min' x max]> with x+ax « 4.1. In this range the numerical solution is very accurate, 
but for larger x + , the last exponential in ( 2.1 3| ) dominates and the numerical solution 



diverges. 

Note that the range [^min' x max] i s sufficient to compare the different approximations and 
the different levels (0,0), (2,4), (2,6), (4,8) around the transition to the non-perturbative 
vacuum, see Fig. || as an example for the tachyon only. Similarly to the level zero solution, 
Fig. |2], the tachyon is initially in the perturbative vacuum <f)(x + ) = 0, then grows exponen- 
tially near x + = 0, slightly overshoots the vacuum expectation value (</>)(2,4)) then settles in 
the non-perturbative vacuum. However using the exponential series ansatz we cannot eval- 
uate the convergence properties around the non-perturbative vacuum. We devote Section 
H| to this issue. 

All three methods agree very well up to the maximum value near x + ~ 1, where the 
non-linear couplings due to the interaction between the various fields become important. 
But even for x + > 1, the solutions are qualitatively very similar. 
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Figure 4: Level (2,4) tachyon condensation: From the perturbative vacuum (j> — we jump to the new 
vacuum. The solution was calculated using the exponential series ansatz with 500 terms and a precision of 
150 digits. 



Let us now consider the solution for the full system of the seven fields at level two that 
we calculated up to x+ ax « 4.1 using the exponential series ansatz (2.13), see Figure |4|. 
The general picture involving the three stages 

perturbative vacuum — > transition — > non — perturbative vacuum 

holds for all fields. The tachyon is drawn with a thick red brush to emphasize its importance 
for the vacuum transition: it is the first component to grow, and thus drives the others 
out of the perturbative vacuum. This affirms the intuitive notion of the tachyon as the 
unstable mode of the string field, indicating the instability of the supporting D-brane. In 
addition to the evolution of the seven fields, we included the vevs in the non-perturbative 
vacuum. In Figure || we zoom on the tachyon in order to show how it oscillates around 
the vev. Typically one would expect that all fields, except the scalars, have vanishing vev, 
preserving translation invariance in the non-perturbative vacuum. But due to the presence 
of the linear dilaton background V ■ x = V + x~, translation invariance is broken in the 
x + -direction. Hence e.g. (B + ) v+=1 ^ is no upset. If V + = 0, then B + can have no vev. 
How exactly the vevs are reached for large x + is studied in detail in Section ||. 

There is more than one way to determine the vevs. The simplest one is the following: 
Simplify the equations of motion (appendix 0) to allow only constant solutions. This 
results in a set of equations where the vevs (cp) , (B^ ) , (F) appear quadratically, all others 
linearly, hence there exist 3 • 2 + 4 = 10 solutions. With mathematica they are found in 
closed form, as expressions depending on roots of 10th order polynomials. Numerically the 
vevs can be evaluated to arbitrary precision. To pick the right solution from the set of ten, 
it is sufficient to consider the tachyon vev: 
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Figure 5: Level (2,4) tachyon condensation: Zoom on the tachyon <j>(x + ) as it oscillates around the vev. 
The solution was calculated using the exponential series ansatz with 500 terms and a precision of 150 digits. 



• Two solutions can be eliminated as ((f)) G C. 

• Six more can be neglected because (<j>) < 0. From the effective potential, Fig. [l|, we 
know it is unbounded for negative field values, hence there can be no finite negative 
vev. 

• One solution has {(p) = 0. The vevs of the other fields vanish, too. This is the 
unstable, perturbative vacuum solution. 

• The last is the non-perturbative vacuum solution: (0) ~ 0.5416. The resulting vevs 
to four significant digits are summarized in Table ||. 



<</>> 


(F) 


(P) 


(B+~) 


(B-) 


(B+) 


(B~) 


0.5416 


0.8808 


-0.1733 


-0.03670 





-0.05190 






Table 2: Vacuum expectation values at level two: the set of constant solutions of the full equations of 
motion corresponding to the non-perturbative vacuum. 

Another way to determine the vevs is to realize that if we expand the string field in 
the universal basis (i.e. using ghost modes and matter Virasoro modes, but no matter 
oscillators), the vevs would be independent of the dilaton gradient, and therefore equal 
to their values in flat spacetime. Using the expression of the Virasoro operators in terms 
of oscillators and dilaton gradient ( |E.6[ ) we can then deduce the vevs of the fields in the 
non- universal basis. Explicitly, at level two in the universal basis and in Siegel gauge, we 
have three fields t, u and v, and the string field in the non-perturbative vacuum is given 
by 

|*) = (t) d|0) + (u) c_!|0) + (v) L m 2 d|0), (2.14) 
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Figure 6: Level (2,6) tachyon condensation: The solutions were calculated using the exponential series 
ansatz with 300 terms and a precision of 150 digits. We have multiplied B ++ by 20 in order to make it 
visible. 



where the vevs (t), (u) and (v) are well known pq |. Plugging the expression ( |E.6| ) for L™ 2 , 
we have 

|*) = (t) ci\0) + (u) c_i|0) + - (v) aLiaiiCilO) - (v) a± x cr x ci\Qi) + -^=V+ (v) aZ 2 \0). 

2 v 2 

(2.15) 

Comparing with our string field expansion Eq. pi] ), and the definition of F ( gg ), we find 

(0> = <t>, </3> = -<u>, (io = i2V5<«>, <s + -> = --^( w ), = -y + («> . 

(2.16) 

And these give precisely the same values as in Table H providing further evidence that our 
action was correctly calculated. 

A curious result worth mentioning is that the following relations 



24 



:F. 



12V2 



i 7 . 



(2.17) 



which, by Eq. ( 2.16| ), should hold for the expectation values, actually hold for all x + at 
level (2,4). This can be roughly seen on Fig. ||. We have checked these relations numerically 
beyond doubt (they hold with a precision of at least 100 digits), but haven't found any 
simple reason to explain them. They must in fact be "accidental" because they do not 
hold anymore at level (2,6). 

Let us now briefly look at level (2,6). The main difference here is that we cannot set 
B ++ to zero anymore. In particular, the equation of motion for B , which at level (2,4) 
contained only terms proportional to B ++ , now contains in particular a term (B ) 2 . We 
show the exponential series solution for the rolling in Fig. 
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Figure 7: Comparison of tachyon condensation at levels (0,0), (2,4), (2,6) and (4,8). The vevs have been 
normalized to unity. The solution at level zero was calculated with the method of steps, and the solutions 
at higher levels were calculated using the exponential series ansatz with the origin of x + chosen such that 
the first coefficient is one. Accordingly the initial condition = K~ 3 was used for the level zero DDE. 
The curves at levels (2,4) and (2,6) are almost indistinguishable. 



Level four solutions Let us now turn to the system of equations at level (4,8). Based on 
the method of conservation laws we used our mathematica code to compute the equations 
of motion. There are now 50 fields to consider, the equations of motion contain 21400 
terms. We want to investigate to what extent the solution changes compared to level two. 
Again we used the exponential series ansatz and computed the first 50 coefficients for all 50 
fields. This calculation took about 15 h on a fast computer. From those 50 fields we decide 
to focus on the tachyon, discussion of the other 49 fields would be redundant, for they are 
all qualitatively similar. Now we compare the tachyon at level two and four, Figure |7| As 
before the vevs of all fields are extracted from the constant solutions. For the tachyon, 
the vev is only slightly bigger at level 4, the numerical value agrees with the estimate in 
1 46]. We notice that the tachyon overshoots more at level four, both on a relative scale 
and in absolute value. As noted in [p0[| , the full solution (no level truncation) converges 
monotonically. One might have expected to see the overshooting decrease monotonically 
with the level as well, since it is a lot more pronounced at level zero than at level two, but 
that is not the case. For completeness, the relative overshooting of the tachyon is shown 
in Table |3[ 



Level 





(2,4) 


(2,6) 


(4,8) 


OAmax -<</>»/<</>> 


18.7% 


7.1% 


6.7% 


8.2% 



Table 3: Relative overshooting of the maximum vs the vev for the tachyon at the lowest levels; it is not 
monotonically decreasing. 
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3. Oscillations around the non-perturbative vacuum 

In this section, we want to consider the equations of motions when the fields are very close 
to their non-perturbative vacuum expectation values. This analysis can be started in a 
relatively general formalism, independent of the level. We start by putting the components 
of the string field into a column vector 

*(a.+) = (<f>(x + ), B + (x + ), B~(x + ), B+-(x + ), . . .) T . (3.1) 

The length n of this vector will be seven at level (2,4), eight at level (2,6), and fifty at level 
(4,8), the highest level that we consider in this paper. Next we write 

¥(x + ) = * + <5*(z + ), (3.2) 

where \l/o is the non-perturbative vacuum expectation value of and 5&(x + ) is a small 
perturbation. We now write schematically the equations of motion for 5^(x + ) at linearized 
level, i.e. keeping only the terms that are linear in 5*&(x + ). Note that because \l/o is a 
solution of the full equations of motion, there will be obviously no constant term in the 
linearized equations. The general structure of these equations is easy to understand. On 
the left hand side, we will write the contributions from the kinetic term; this involves the 
fields at light-cone time x + (not retarded) and at most first derivatives of the fields. On 
the right-hand side we will write the contributions from the interaction term. These all 
involve the same retarded light-cone time x + — 7, and at most 3L derivatives, where L is 
the level. This is so because the number of derivatives is at most equal to the total number 
of Lorentz indices carried by the interacting fields, and each field of level L can carry at 
most L indices. So the general form of the equations of motion is 

8&(x + ) + A5&(x + ) = 

B <5*(x+ - 7) + B x 5V'(x + - 7) + B 2 <5*"(x+ - 7) + . . . + B 3L <5* (3L) (x + - 7), (3.3) 

where all the details are hidden in the n by n matrices A and B m , m = 0, . . . 3L. We can 
now make the ansatz 

5V(x + ) = e ux+ 3, (3.4) 

where uj is a complex number and S is a vector of complex numbers. Since the equations 
Q3.3D are real and linear, we can always find a real solution by adding its complex conjugate 
to fl3.4p . Note that with this ansatz, all fields oscillate with the same frequency Imw with 
exponentially decaying or growing amplitudes, according to the sign of Rew. The relative 
amplitudes and phase shifts between the fields are encoded in 3. Plugging this ansatz into 
Eq. ( |3.3| ) and multiplying by e 7 ^, we obtain the equation 

(we 7U / + e 7 " A - B - u B 1 - . . . - uj 3L B 3L ) 3 = 0, (3.5) 

where I is the n by n identity matrix. After naming 

M(w) = we 7W / + e 7 " A - B - uj B 1 - . . . - uj 3L B 3L , (3.6) 
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Figure 8: The distribution diagram of det M(uj) at level (2,4). It has three segments, L\ (retarded), L2 
(neutral), and L3 (advanced) with respective slopes —, 0, and ——. 

we see that Eq. fl3.5|) has a nontrivial solution for S if and only if uj is such that 

detM(w) = 0. (3.7) 

It is easy to see from (|3.6| ), that det M{oj) is an exponential polynomial in uj, i.e. a function 
of the form 

n 

det M(oj) = J>» eft " Q = po<fa.<...<Pn, (3.8) 

3=0 

where Pj(uj) are polynomials of degree dj. In our particular case, we have /3j = j'7. 

So we are interested in finding the roots of an exponential polynomial (|3.7|), and partic- 
ularly in the signs of their real parts because they will determine whether the perturbation 
5^(x + ) will eventually die out or not. This problem has been studied quite extensively 
(see for example |^3|] ) . The general answer is that one can write expressions for the asymp- 
totic values of the zeros (i.e. the zeros u with \u\ —> 00). More simply, one can use the 
distribution diagram to approximately locate the zeros. Namely, we plot the points Pj with 
coordinates (f3j,dj); we then draw the convex polygonal line L that joins Pq to P n , such 
that its vertices are points of the set Pj and such that no point Pj lies above it. This is 
illustrated in Fig. |8] for the case of level (2,4). Now let us denote the successive segments 
of L (from left to right) by Li, L2, . . . L/-, and let m,^2, • • • Mfc denote their slopes. Now 
the general result is that there exist positive numbers c\ and C2 such that all the zeros 
with norm greater than C2, are located inside the union of the strips V r defined by 

|Re(w + fj, r logw)| < c\. (3.9) 

It is clear that, for large |cl>|, these strips are located in the left half-plane if \i r is positive 
and in the right-half plane if \x r is negative. If /i r = 0, the strip is vertical and contains the 
imaginary axis. At level (2,4), we have three strips with \i r respectively equal to ~, 0, and 
— -. These strips are sketched in Fig. ^ £). One can say a little more. In particular, in any 
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Figure 9: Location of the zeros of detM at level (2,4). i): The asymptotic analysis tells us that there 
exists a C2 > such that outside the circle of radius C2, all the zeros are located inside the retarded strip 
Vi, the neutral strip V2 and the advanced strip V3. ii): Roots found numerically for —200 < lmoj < 200. 
The blue dots are simple roots while the red dots are double roots. 



region R denned by 

I Re(w + fj, r logw)| < c\, I Im(w + [i r log u) — a\ < b, (3.10) 
with no zero on the boundary, the number n(R) of zeros in R satisfies 

l-n r + -P<n(R)<-P + n r -l, (3.11) 

7T 7T 

where n r is the number of points of the distribution diagram on L r and (5 is the difference 



in values of j3j between the end-points of L r . In Eq. (3.10), a is an arbitrary real number, 
and b is an arbitrary real positive number. So, for example with [i T = (neutral strip), 
Eq. ( j3.10|) defines a rectangular box of width 2ct and height 2b centered on a X i. Therefore, 



Eq. (|3 -11 ) just means that the average vertical density of the zeros along the neutral strip 



is given by j3/ir. For the other branches with /x r 7^ 0, the picture is almost the same, 
just a little twisted. But in the limit |u>| — > 00, log(cj) is going to be roughly constant 
in the box, so ( |3.10 ) defines an almost rectangular box centered on —fi r \og(u) + a x i. 



This means in particular that each strip contains infinitely many zeros with arbitrarily 
large norm. And it has the following immediate consequence: At level (2,4), detM(cj) 
has infinitely many zeros with arbitrarily large positive real part. This means that if we 
specify a generic initial condition for the string field ^(x + ), < x + < 7, with *&(x + ) close 
to the non-perturbative vacuum, then we should expect that it will start to oscillate with 
exponentially growing amplitude, at least until the linear approximation becomes invalid. 
For concreteness we show, in Figure I ii), the zeros of det M(uS) found numerically in the 
range —200 < Imw < 200. The picture agrees with the asymptotic analysis. We see in 
particular that the zeros are aligned along branches in each of the strips, and that along 
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Figure 10: Distribution diagrams of some simple linear delay differential equations of i) retarded type, 
ii) neutral type, and Hi) advanced type. 



these branches the spacing between roots is asymptotically constant. Interestingly, in the 
retarded strip there is a branch of double roots (the double zeros are plotted in red in Figure 
|9| ii)). It turns out that this is related to the fact, already mentioned in Section ||, that at 
level (2,4) there are only five independent fields because of the relations ( 2.17| ). 



In Fig. g we have denoted the segment with positive slope as retarded, the one with 
zero slope as neutral and the one with negative slope as advanced. This denomination is 
easily understood by looking at the following three simple linear delay differential equations 
(we refer the reader to Appendix [A] for definitions and more details on delay differential 
equations). First, suppose we have an equation of the form 

f'(t) = f(t-l). (3.12) 

This is an equation of the retarded type since the right-hand-side includes / at the retarded 
time t — 1. Making the ansatz f(t) = e ut , we obtain the characteristic equation for u 

toe" -1 = 0, (3.13) 



whose distribution diagram is shown in Fig. 10 Q. It has one segment of positive unit slope; 



this justifies the denomination "retarded" for such segments. Now we consider the neutral 
delay differential equation 

f'(t) = f'(t-l), (3.14) 

which has for characteristic equation ue w — oj = 0. Its distribution diagram, shown in 
Fig. |l0|ji), has one segment of zero slope, hence the denomination "neutral" for such seg- 
ments. At last, the equation 

f'(t) = f(t + l) (3.15) 

is of the advanced type because the right-hand side involves / at the advanced time t + 1. 
Its characteristic equation is u — e u = 0, whose branch diagram is shown on Fig. |l(]m). It 
has one segment of negative slope, which we therefore name advanced. 

To summarize, at level zero the linearized equations of motion are purely of the retarded 
type because the characteristic equation for the tachyon is [fjOf: 



(w - l)e 7 ^ + 2 = 0, (3.16) 



and its distribution diagram is thus as shown in Fig. Ki). From the above discussion, we 



therefore know that its roots of large absolute value all have negative real part, and it could 
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thus have at most finitely many roots with positive real part. It turns out that all the roots 
have negative real parts, and the motion around the non-perturbative vacuum is stable at 
this level. However, we have shown that at level (2,4) we obtain two new segments, a neutral 
one and an advanced one which ruins the stability around the non-perturbative vacuum. 



It is interesting to ask what happens at higher levels. We show, in Fig. 11, the distribution 
diagram found at level (2,6). Note that at this level, it is not anymore consistent to set 



d; 




Figure 11: The distribution diagram of detM(o;) at level (2,6). It has four segments L\, . . . ,La with 
respective slopes -, 0, —-, and — -. 



B ++ (x + ) to zero, so we have a total of eight fields. One of the most striking differences 
with level (2,4) is that there is one more advanced segment. This can be understood from 
the fact that the action at level (2,6) contains higher derivatives at retarded times. One 
notices also that the retarded segment L\ has become much shorter. With our code, we 
were able to calculate the action and equations of motion up to level (4,8). From this large 



calculation we show in Fig. 12 the result for the distribution diagram. At this level we have 




Figure 12: The distribution diagram of detM(oj) at level (4,1 
respective slopes -, 0, — — , — -, and — -. 



It has five segments Li,...,Ls with 
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fifty fields. It turns out that the calculation of det M(uS) is a difficult problem. With the 
definition z = e 7a; , we see that the entries of the matrix M{uS) are polynomials in the two 
variables u and z. We explain, in Appendix why the computation of the determinant 
of a polynomial matrix is hard and how to overcome this difficulty. 

Interestingly, we see that at level (4,8), we have one retarded branch, one neutral 
branch, and three advanced branches. In this sense, the equations of motion become more 
and more advanced as the level is increased. We can actually express this idea differently; 
one could ask if we can write an effective equation of motion for one field, e.g. the tachyon 
(f)(x + ), that would give rise to the same characteristic equation as Eq. Q3.7| ). The answer 
is very easy; indeed if we make the ansatz (f)(x + ) = e wx+ , it is immediately clear that the 
equation of motion 

n 

Y,Pj(d + )4>(x + +jj) = (3-17) 

is precisely the same as Eq. ( ^!?| ) , where the pj 's are the polynomials defined by Eq. (^) . 
We now do the usual splitting of this differential equation by putting one of the terms 
with highest derivatives on the left-hand side and all the other ones on the right-hand side. 
So we choose a j max such that the degree dj max of the polynomial p Jmax is maximum (i.e. 
4?' m ax ^ dj, j = 0, . . . n). If the distribution diagram has a neutral segment, then we can 
choose, for j max , any j along this segment. After subtracting 7 max 7 from x + , the equation 
of motion thus takes the form 

P; max (3+) </>(*+) = - E P;(d+)<Xx + + (j-j max ) 7 ). (3-18) 

j^j 'max 

What we can immediately read from this equation is that, at level higher than zero, the 
non-locality is not anymore concentrated at one point in the past x + — 7, but instead takes 
all the values {j — j max ) 7, J ' = 0, . . . n. While this precise range depends on the value of 
imax that we chose, it is clear that it will extend in the future by, at least, an amount 
^future 7 equal to the length of the projection of the advanced segments on the /3j axis. 
And similarly, it will extend in the past by, at least, an amount n pas t 7 equal to the length 
of the projection of the retarded segments on the (3j axis. 

If we specify generic initial data for DDEs (|3.18| ) whose lowest derivative appears at an 
advanced time (which we call advanced DDEs), we should not expect that smoothing will 
occur when solving with the method of steps. In fact, we show in Appendix |A| an example 
of an advanced DDE which does not even possess a continuous solution for a given initial 
data. For a system of advanced DDEs, on the other hand, smoothing may occur. For our 
rolling solution, we require an analytic solution. This is not in conflict with the preceding 
remark because we are not free to choose any initial data, as we have shown that the initial 
conditions are essentially uniquely determined by demanding that the tachyon lives in the 
perturbative vacuum in the infinite light-cone past. 

It is clear from the distribution diagrams in Figs. ^ and [l^, that nf uture and n past both 
increase from level (2,4) to level (4,8). Going from level (2,4) to level (2,6), on the other 
hand, we see that although nf uture increases, n past decreases. But if we look only at levels 
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(L, 2L), we nevertheless still expect that both nf uture and n pas t will continue to grow beyond 
level (4,8), and that in the limit L —> oo, the non-locality extends over the whole discrete 
range {re 7, n£Z), 

At last, we remark that at infinite level, if we subsequently take the limit V + —> 
0, which is equivalent to shrinking the delay 7—7-0, the range of non- locality becomes 
continuous and extends over all light-cone times. We thus recover, in this limit, the non- 
locality characteristic of usual (zero dilaton) string field theory. 



4. Discussion and Conclusions 

We have identified the equations of motion for the open string tachyon rolling in a linear 
dilaton background as being delay differential equations. This class of equations has been 
widely studied in the mathematics literature. Although no method to analytically solve 
non-trivial DDEs is known, the method of steps is in general an efficient numerical method 
(and essentially the only method) once the initial data is provided. We have used this 
method in order to solve the equations of motion at level zero, and found that it is both 
efficient and accurate; in particular our solution agrees very well with the results from an 
exponential series ansatz, and can be continued much further in light-cone time (Fig. ^). 

At truncation level two, we found that the method of steps can be used only when 
we further truncate the action by either removing all higher derivatives, or by keeping 
only the scalar fields. But we saw that the results from these approximations all agree 
reasonably well with the more accurate solution from the exponential series ansatz, at 
least until the first maximum has been reached (Fig. ^). The obtained numerical solutions 
can be continued for much larger light-cone time, where it is seen that they converge 
towards the non-perturbative vacuum. If we don't simplify the equations of motion, our 
numerical solver gives us wrong solutions. We have understood that this problem is due 
to the fact that our equations of motion possess derivatives of order higher than one at 
retarded times. In general, such equations are not even guaranteed to possess continuous 
solutions. While we still expect the equations of motion to have an analytic solution for 
the initial data corresponding to a tachyon in the perturbative vacuum in the far light-cone 
past, the method of steps is rendered unstable by the higher derivatives. 



We have already mentioned that the diffusion equation method [36, 37, 38] has been 
successfully used to solve numerically the light-like rolling problem at level zero ]0q| . It 
is then natural to ask wether this method can give more stable solutions at higher levels. 
The diffusion method has the virtue of applying to a larger class of non-local equations. 
But in the case of DDEs, we believe that it is completely identical to the method of steps. 
It would then face exactly the same problem as the method of steps at higher levels. To 
illustrate this, we briefly review the diffusion method applied to the equation 

<f>'(x + ) -(f>{x + ) = -cf)(x + - j) 2 . (4.1) 

First we define a function of two variables x + and r by <3?(x + , r) = e ir9+ <j){x + ) = <p(x + +rj). 
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This definition is implemented by the simple diffusion equation 



<9+3>(x + ,r) = -d r <5>(x + ,r). (4.2) 

7 

The idea is then to solve the diffusion equation on the region x + > and < r < 1. 
The subtle part is to determine the boundary conditions. It turns out that the boundary 
conditions on the r = 1 axis are determined by the equation (which is now local in x + ): 

$(x+, 1) - ~ [d r $(x+,r)) r=1 = $(x+ 0) 2 . (4.3) 

A simple algorithmic resolution of the diffusion equation is the following: 

1. Specify boundary conditions on the segment x + = 0. This is the initial data. 

2. Solve the diffusion equation ( |4.2| ) with the given piece of boundary conditions. One 
can easily see that this amounts to copying "diagonally" the values on the segment 
x + = to the horizontal segment defined by r = and < x + < 7, see Fig. [D^. 

3. To solve the diffusion equation further, one needs first to specify the initial conditions 
on the segment defined by r = 1 and < x + < 7 using Eq. (|4,3|). Since the simple 
diffusion equation has been solved exactly in the previous step, replace ^ by d+, 
and the right-hand side is simply the initial data. One obtains precisely the original 
equation Q4.1| ) (with 4>(x + ) replaced by <l?(x + , 1)). 

4. Solve again the diffusion equation by copying the obtained boundary values on the 
segment r = 1 and < x + < 7 to the segment r = and 7 < x + < 27. 



5. Solve again Eq. (|4.3|) on the segment r = 1 and 7 < x + < 27. From the same 



argument as above, plus the fact that one can replace the right-hand side of ( [4.3|) by 
<£(x + — 7, l) 2 , this is solving the DDE equation ( fl.l| ) on the interval 7 < x + < 27 
given the values of (j> on the interval < x + < 7. This is precisely the method of 
steps as described in Appendix |A} 

The algorithm then continues by repeating points 4 and 5, which is equivalent to solving 
the original DDE ([O]) by the method of steps. 

To summarize, the diffusion method applied to a DDE "mimics" the method of steps. 
We cautiously note, however, that a different numerical method to solve the diffusion 



equation may change the above claim if Eq. (4.2) is not solved exactly as in our algorithm. 



In order to find accurate numerical solutions taking all derivatives into account, we had 
to resort to the exponential series ansatz. This method generalizes obviously from level 
zero; in particular one still can calculate recursively the coefficients of the series. The 
problem with this method is that the number of coefficients needed grows exponentially 
with the light-cone time x^, ax until which we want our solution to be accurate. Moreover, 
we must calculate the coefficients with many digits precision in order to account for the 
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Figure 13: Step 2 in the algorithm to solve a DDE in the diffusion equation approach. The initial data 
at x + = is copied "diagonally" to the r — boundary. 



numerical cancellation errors. At levels (2,4) (Figs. |3] and |5|) and (2,6) (Fig. |6|) we were 
nevertheless able to find solutions that are accurate up to a point where the string field is 
convincingly seen to converge towards the non-perturbative vacuum. Moreover, we have 
found that the convergence (as measured by the overshooting of the first maximum) is 
substantially faster at levels two and four than at level zero. At level four, it is a little 
bit slower that at level two (Fig. |?]). In order to compare this to the analytic solution of 
Hellerman and Schnabl [|30]], it is instructive to look at the tachyon component of their 
solution (in Schnabl gauge). Normalizing the tachyon vev in the non-perturbative vacuum 
to one, it is 

A ( °V) = 77-^' ( 4 - 4 ) 

1 + e x 

a monotonic solution. So the fact that our solution at level two has become more monotonic 
is a good indication that level truncation yields qualitatively the same solution (in Siegel 



gauge). We note, however, that if we expand (4.4) as a power series in e x (what we have 
called an exponential series), it has a finite radius of convergence; the series diverges for 
x + > 0. On the other hand, the exponential series that we obtained have a much larger 
radius of convergence, perhaps infinite as can be shown at level zero from the asymptotic 
behavior of the coefficients of the series [|3C(| . It is worth mentioning a curious fact. Had we 
taken the level-zero equation of motion (with vev normalized to unity) cj)'(x + ) — 4>(x + ) = 
—(f)(x + — 7) 2 , and dropped the delay 7, we would get 

<p'(x + ) -<p(x + ) = -(p{x + ) 2 , 

which has the solution 

0(x +x 



1 + e a 



Precisely the same as (4.4)! This might just be a coincidence; it could also point at 



something deeper which may be worth examining in a further work. 
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To sum up, our numerical solutions at levels two and four strongly indicate that the 
string field will eventually converge to the non-perturbative vacuum, in agreement with 
the analytic solution of |3(|. On the other hand, we have also shown that, at levels higher 
than zero, the equations of motion, linearized around the non-perturbative vacuum, admit 
infinitely many oscillating solutions with exponentially growing amplitude. This suggests 
that the string field should not converge to this vacuum. These two apparently contradic- 
tory results can be reconciled in different ways. Below we describe two possibilities. 

The first possibility is that the string field somehow avoids the infinitely many expo- 
nentially growing modes. This is certainly possible because the initial conditions (i.e. the 
string field sitting in the perturbative vacuum in the infinite light-cone past) are very spe- 
cial: With this point of view, we should conclude that the motion is unstable: should an 
external perturbation, or a quantum effect, move the string field the tiniest bit away from 
this particular solution it would not converge. We should also stress that our analysis of 
the linearized equations of motion doesn't tell us how non-convergent the rolling would 
be. It might approach a limit cycle (in phase space) as is sometimes the case when the 
linearized equations of motion have a finite number of growing solutions. Our situation, 
however, is different because we have infinitely many growing solutions in the vacuum. 
This could render the rolling diverging. Or it might not be that bad; after all, we could 
imagine expressing a well-behaved converging function as a series in the infinitely many 
diverging modes ||43| ]; this would not be possible in the finite case. We emphasize here 
again the special character of the initial condition of the tachyon sitting on the top of the 
hill in the infinite light-cone past, which we have shown uniquely fixes the initial data up 
to a shift in light-cone time. The method of exponential series, which gives a convergent 
solution, can work only with this particular initial condition, which is also the initial con- 



dition of the analytic solution of [30]. For a generic initial condition the coefficients cannot 
be determined recursively, because "negative" modes have to be taken into account. In the 
level zero case, with (j)(x + ) = Yln°=-oo a n£ nx+ , one obtains 



K 3 



-In °° 



a„ 



n 



— a m a n - m . (4.5) 



In other words: in order to calculate one coefficient, we need to know all coefficients, 
including the one we are interested in! The same problem is found in the DDE approach: 
to find initial data that are consistent with the equation of motion, we need to solve the 
very same equation of motion. This is a simple problem only when we require that the 
tachyon be in the non-perturbative vacuum in the infinite past. 

We leave it for further work to study different initial conditions. At level zero it was 
done in [35|; but at higher levels we expect different results due to the higher derivatives 



and divergent modes. 

A second possibility to reconcile a convergent rolling solution with infinitely many grow- 
ing modes is that these modes do not satisfy the "out-of-gauge" equations of motion 6 . 



°We thank T. Erler and M. Schnabl for emphasizing this point to us. 
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Indeed we have derived the equations of motion from the gauge-fixed action. This means 
that for a solution we are guaranteed that (<J>, Q^?o + * *o) = for all $ in Siegel 
gauge. But it is possible that Q^o + itself doesn't vanish. That such solutions 



can exist was shown in 47 



The two different explanations above differ in particular by the interpretation of the 
growing modes. In the first possibility, they do have some physical meaning and make the 
rolling unstable. While we do not expect to find open string degrees of freedom in the 
non-perturbative vacuum, it is not unreasonable to think that the unstable modes could 
correspond to the closed tachyon instability. After all, we have coupled the open string 
field theory to the closed sector, although only via the dilaton. In the second possibility 
mentioned above, however, the growing modes are simply unphysical because they do not 
obey the out-of-gauge equations of motion. In order to find out if this is the case, we would 
have to derive the equations of motion from the gauge-unfixed action and verify whether 
the growing modes are still solutions of these full equations of motion. We hope to be able 
to report on this in a future publication. The present paper gives only a partial answer to 
the question of the physical meaning of the growing modes. At the end of Section ||, we 
wrote an effective equation of motion for one scalar field, whose solutions are precisely the 
exponential modes solutions of the whole system of equations of motion linearized around 
the non-perturbative vacuum. We have then found that the decaying modes correspond to 
positive delays in the effective equation of motion, thus non-locality in the past, and that 
growing modes correspond to negative delays, i.e. non-locality extending into the future. 
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A. Delay differential equations 

Differential equations describing real life systems like chemical reactions, predator-prey 
relations, etc... often need to incorporate the state of the system not only now, but also at 
times in the past. A simple example is the following delay differential equation (DDE) 

y '(t) = -y(t-l) (A.l) 

with the constant delay term y (t — 1). Notice the non- locality: the derivative is needed at 
t, but the function itself is evaluated a finite distance away at t — 1. This equation is of the 
retarded type. In order to make this a well posed problem it is not sufficient to look for 
a solution starting at a point to with the value of the function y (to) = yo fixed as in the 
Cauchy initial value problem for an ordinary differential equation (ODE). Instead one has 
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to supply the initial data <p (t) over a continuous range of points. Choosing for example 
4>(t) = 1 one has as a well posed problem 

(y'(t) = -y(t-l) t>0 

\y(t) = l t<0. 

In contrast to an ODE problem the solution to the above problem has discontinuity points 
of order n where the n-th derivative is discontinuous. One can quickly see how this comes 
about. Approaching from the left, lim^g- y' (t) = 0; however approaching from the right 
one finds 

lim y '(0) = -y(-l) = -l. 
t->0+ 

This discontinuity propagates further along the solution. The second derivative at t = 1 
reads 

y" (1) = -y' (0) , 

and is thus discontinuous. However, the first derivative at t = 1 is y' (1) = — y (0), which 
is continuous. Hence the solution is of class C at t = 1. This process continues at every 
positive integer t = n, where the solution is C n ; hence the discontinuities are smoothed out 
as t increases. The propagation and existence of discontinuity points is a generic feature of 
DDEs. In general the solution only links continuously with the initial data at the starting 
point to- 

Consider as a second example a DDE of the advanced type 

(y'(t) = -y(t+l) t>0 

W) = l t<0 . 

The only variation to (|A.1| ) is that the derivative is taken at the present instant, but the 
function is needed at a future point. For — 1 < t < we have y (i + 1) = y' (t) = and 
similarly the solution vanishes for all positive t. The solutions to both examples, ( |A.1| ) 
and ( |A.2| ) , are plotted in Figure [l4|. The discontinuity of the advanced solution at t = is 
noteworthy: for this type of equation only special initial conditions allow for a continuous 
solution. In our example that special choice among the class of constant initial data is 
y(t) = 0, t<0. 

The above examples can be generalized further to include any number of delays, pos- 
sibly variable or time dependent delays r (t) or even state dependent delays T(t,y(t)). If 
there are derivative terms at delayed points, y' (t — a (t, y (t))) the equation is said to be of 
the neutral type ('NDDE). In general the delay functions T(t,y(t)) and a(t,y(t)) can be 
completely independent. If / depends on y over a continuous range, the problem is called 
a Volterra functional differential equation, but for the present purpose, one can restrict to 
only a finite number of points in the past required to compute the current rate of change 
in the system. 

The general definition of a first-order delay differential equation problem then reads 

' y'(t) = f (t, y(t),y(t - n(t, y(t))), ...y(t- r m (t, y(t))), 

y'(t-a 1 (t,y(t))),...y'(t-a n (t,y(t)))) t <t<t f 
y(t)= <f>(t) t<t . 



-30- 



2 




_2 Li i i i i i i i i i i i i i i i i i i i i i i_ 

2 4 6 8 10 

t 

Figure 14: Solutions to the retarded and advanced examples with identical initial data. The retarded 
equation possesses a unique continuous solution for any initial data, but the advanced solution is in general 
discontinuous. 



Extending to systems of DDEs or higher order derivative terms is then straightforward. 
There are further major differences between ODEs and DDEs in addition to the above 
discussed discontinuities. 

• Different initial data <ft(t < to) can give rise to the same solution y (t > to) 

• The solution may be non-unique in the state dependent case 

• The solution may terminate at a point 

• In the neutral case generally no smoothing of discontinuities occurs 

• Bounded solutions may oscillate or show chaotic behavior even in the one-dimensional 
case, while ODEs can oscillate only in at least two dimensions and behave chaotically 
in at least three dimensions (Poincare-Bendixson theorem). In fact oscillations are 
typical of DDEs 

• infinitesimal delays introduced into ODEs can stabilize or destabilize the solution 

• For scalar NDDEs no smoothing occurs, but for systems of NDDEs smoothing may 
take place 

Existence and Uniqueness Theorems on local existence and uniqueness have been 
developed for very general assumptions on the delay functions r(t,y(t)) and a(t,y(t)). 
But for the case at hand in this paper one can use familiar results from ODE theory. 
Assume there is only one delay, a positive number on the interval of interest [tojfy] 

t = a > 0. 
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This case of constant delay then reduces to the following ordinary differential equation ini- 
tial value problem on the interval [to, to + T ] between the initial point and the corresponding 
point shifted by one delay 

f y' (t) = f (t, y (t) , (t - t) , (/)' (t - r)) t < t < t + t 
\y(t ) = <p(t ). 

Standard results then guarantee local existence and uniqueness of a continuous solution in 
[to, to + 8] for some 8 > provided that / (t, y (t) , 4> (t — r) , eft' (t — r)) is continuous with 
respect to t, Lipschitz continuous with respect to y, <f>, 4>' and that 4> (t) , <j>' (t) are themselves 
Lipschitz continuous. Note however that there exists a whole family of discontinuous 
solutions, linked to the unique continuous one by a constant offset. 

Method of steps In order to obtain the solution on the full interval it is convenient both 
theoretically and computationally to iterate along the intervals [to + j ■ r, to + (j + l)r], 
j = 0,1,2,... using the solution calculated for the previous interval. This is known as 
the method of steps, the principal algorithm to solve DDEs numerically. It is applicable 
in the case where all discontinuity points are known in advance, e.g. for a constant delay. 
Solving the DDE then reduces to solving ODEs on the intervals between two discontinuity 
points using standard ODE methods such as linear multistep solvers and gluing the solution 
together at the interval boundaries. For continuing over several such intervals the ODE 
method has to supply a continuous (polynomial) interpolation of the solution which matches 
the accuracy order of the solver. For concreteness, let yo be the solution on the interval 
[to, to + t], then for the solution yiin the next step the ODE reads 

U (t) = / (t,yi (t) ,?/o (t-T),y' {t-T)) t + r < t < t + 2r 
[yi (to + r) = y (to + r). 



Further references are, for example, [J43|, [44] and [45]. 
B. Lagrangian at level (2,4) 

We present the full Lagrangian density with fields up to level two and interaction terms up 
to level four in the linear dilaton background with light-like dilaton gradient V 2 = 0. The 
Lagrangian can be split into two: one is explicitly ^-dependent and vanishes for V = 0, 
the other survives for V = 0. While V could always be eliminated using the modified 
momentum conservation ^ k{f + iV^ = 0, we keep it explicit for easier comparison with 



1 41]. The interaction term can be split up according to the level n and the number of 



derivatives k, C^ 1, Overall we have 



S = -^Jd D x e- v -*^-{£ Km + C v Kin ) + l -K* (£ Int + CYnt)] 
C M = C^ + £(2) + + L f + C f + Cf + Cf. 
We use the following expansion of the string field into Fock space states 
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Working in the light-cone, we have 



and for example 



e -V-x = e V+x~ 



= -V+B-. 



We split the Lorentz indices into light-like and ordinary components 

M=(0,l,2,...D-l)->. (+,-,2,3,... D-l) = (+,-,»). 

In order to simplify the equations of motion, we need not consider all 379 independent 
components of B^ U ,B U , (3,(f) from the expansion (where B^ v is symmetric). Several facts 
aid in reducing the number of fields: 

1. Any symmetric tensor can be reduced to a traceless and a pure trace part: Bij = 
Sij + ^Frjij with Sij traceless. 

2. For any field <3?i € ft = {B\ B +t , B lj , B~\ S tj }, it only appears in the action 
coupled to another field $2 £ Si> Thus we can safely set these fields to zero and 
trivially satisfy their equations of motion, but still allow for a non-trivial solution for 
the remaining fields. 

We end up with the following eight fields to solve for: 

{</>, B+, B-,B ++ , B + -,B-,F, (3). 

Fields with a tilde are defined by 

Explicitly the kinetic and interaction C without the extra dilaton contribution contain the 
following 41 terms: 7 



T ~For brevity we set a' = 1. 
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c,Hn = d^4> + d^B v d»B u + d^B vX d»B uX - 

-<j? + B^B" + B^B^ - f3 2 
= ^B^-^V-\^d^- ^dMB^ + ^ V ^ ^^ t ^ 

4 fl = ^W^-^M^^M^/ + ^^ 



4 4) = ^ * p ' + 5 ^ / - ^4>B ^B p o 



mo ^19 ^19 

- 243 V B ^ B MA + 243 W £ A + — „ A 0"i? ^ A 

" ^ ^ + ^ 5 W + *|^<WW B 

^d^B^d^ 

32 1 28 

—d^dpd^ B pc B t» - —0 d^dpB ^d a B >"> 

1 28 1 28 

- m ?^d v B'"d p d IT B ^ + —d^d p B a »d*B *> 

In addition we obtain eight interaction terms and four kinetic terms from introducing 
the dilaton. The kinetic terms however do not contribute to the equations of motion, they 
are essentially total derivatives. 



C v Kin = V v B»d v B p + V v B^d v B^ + v v 4>d v ^ - v u pd v p 

-f^d v B^ - ^V,B»B/4> ~ W 



Writing the action explicitly in terms of these eight fields gives the Lagrangians: 



£(2) 



9\/2 



-2B + ~ + F 



) 



11 



2^ 
3 



d+B+ + d-B- 



16 

9y/2 



a 



dl4>B ++ + 2d+d-4>B + - + d 2 4>B—^j 
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81 ^ 486 V J y 

+ ™j>( 2B++B— + 2B+-B+- + -F 2 ) + (-25+~ + P) + 

243^V 24 / 243 HY \ J 8V 7 
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+B-~d 2 (-25+~ +f)} 

+ (-25+- + f) (d 2 <^5++ + 2d + d^B+~ + d 2 jB~-) 

- ^c/ jd+05++<9+ (-25+- + f) + d-4>B+-d + (-25+- + f) 
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+ ^"'9+^ (-B— a_5++ + B+-8-B+- + B+-&+B++ + 5++d+5+-) 
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- |^c/ (9 + <^5++5+- + d + d-4>B ++ B— + d + d^B + -B + - + d 2 <^5+- 5") 

- (5+5— a_5++ + d+B+-d-B + - + d+5++d+5+- + 
88- \/2 , 



243 



c/ (d+^5++d + /3 + d+05+~d_/3 + d_05+-<9 +/ 3 + d-^B~d-Pj 



+ (5++<9 2 /3 + 2B + ~d + d-p + B~ 9 2 
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- ||a /2 \d + 4>d + B ++ d\ B++ + d-^d-B" d 2 B— + d + 4>d+B— d 2 + £++} 
_ l^' 2 |5_^_S++a 2 B~ + 3+<^ (d-B ++ d\ b+- + d-B—d 2 _B+-^ } 

- ||a /2 {<9_0 (d+i? ++ d 2 B+- + d+B— d 2 B + -^j + 2d+4>d + B + - d+d-B ++ } 

- ^a' 2 ^2d + 4>d-B + -d + d-B + - + 2d-4>d + B + - d + d-B + ~ 

+2d-j>d-B + -d + d-B—} 
+ ^a' 2 j a 2 (a +J B++) 2 + a 2 (a_B~) 2 + d 2 <£ (5_B+-) 2 + d 2 J (d + B + -) 2 } 
+ |^a' 2 a 2 ^ + 5-a_5+- + |§a' 2 {2d 2 iM+-3_5++ 

+ 2<9 + <9„0 (d + B + ~d + B ++ + d-B + ~d-B— + d+B—d-B ++ 
+d+B + -d_B + - 
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11 

¥ 

121 , 

243 



242 



243 .a>V+B-(-2B+- + P)^ + ^V+B-M 



_ \ 2 _ 44 _ _ / _ _ \ 

+ — o' ( i ' ' b J + — n'r • b u (<r /r + .9 b J 



! ' (>V ^a /3/ V + (d + B-B ++ d + ^ + d + B-B + -d-4> + d-B-B + -d + 4> 



243 

+d-B~B—d-fj 

88\/2_ / 3/2 T ^+ / .,2 ,, P 



+ ^f«' 3 /2 F+ (dlB-B ++ 4> + 2d+d-B-B+-~4> + &_B-B—$) 



88\/2 . 

H —a' 

243 



3 / 2 V+B _ (b++<9 2 <?!> + 2B + ~d + d-4) + B—d 2 ^j 



c Kin = - d + <t>d-(j) - (d+B ++ d-B— + a+s— d_B ++ + 

- {d+B+d-B- +d+B-d-B+) + d+/3d_/3 

+ ^7 |-0 2 + 2B++B— + 2 + ^ - 2B+B~ - /3 2 j 



C. Equations of motion at level (2,4) 

With the solution B ++ = and the retarded point y + = x + — 2a'V + log (K), the seven 
remaining equations of motion become: 



Equation of motion for the tachyon tfi (48 terms): 

= - 

a' 

( cc en in 

+K 3 |^v / 2>( y + )/3(y + ) " + ^V2^F{y+)B+\y+) 

- A^F(y+)0(y+) + 139 f 3 ( 7 f )2 - ^v / 2 / 3( y +)B+-( y +) 
+ ^V^/3(y+)B+ V) - |/3( y +)^( y + ) + ^/3(y + ) 2 

- |jjvW^B+( y +)B+-V) - g^v / 2v / a^+-( y +)B+V) 
+ H^2i?+-(y+)0(y+) + ^B + -{y+) 2 - ^cf>(y + )B + ' (y + ) 

+ ^*>B+'{y+f - |^5 + ( y + )5-(y + ) + 0(y + ) 2 } 
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+K 3 V+{^a>B+-(y+)F>(y+) + ^ a 'F(y+)B+- , (y+) + ^V2V^V)£ 
+ ^V2a'B+-(y+)(3>(y+) + V2a'/3(y+)B+- '(y+) 



"(v + ) 



729 
242 



+ |g^/3(y+)B-(y+) + ^V2V^V- V)# + V) 
- -^a>(y+)B+-'(y+) - |W-(y+)S+- V) 



27 

+ ^^B+-( V +)B+\y+] 



122 
243 



^^ J B+-(y+) J B-(y+; 



Q9 256 
- - —V2V^B-(y+)B+(y + ) 

44 22 
+ -=- a 'S-(y + )5 +/ (y + ) - — V 



243 
40 



*'£r(y+)0(y+)} 
+K\V+f{^ a 'F{y+)B-~(y+) + i| v^i9(y + )B-(y+) 



+ ^'^-(y+^-'V) + ^V2V^ 3 J B-(y+)^-'(y+) 
+ ^|^^ 3 B+-(y+)i?- / (y+) + ^V2V^ 3 B-(y+)B+'(y 



176 



J B-(y+)B+- , (y + ) + a> 2 B + -(y + )B-'(y + ) 



- |V2a' J B--(y+)0(y+) + |^'£-(y+) 2 } 

128 ,9 _ . /. , . 256 /2o+ _ 

729 

+ ^^^"(y+)i?-(y+)} 
+^ 3 (^ + ) 4 { + ^« /2 ^-( y + ) 2 } 



Equation of motion for B ++ (30 terms): 



0= (x + ) + 

40 
729 



a' 



+K 3 {^a><j>(y+)F"(y+) - ^a> 'F' '(y+)<f>> '(y+) + ^«'*V)W 



+ ^^a>(y + )/3"(y + ) - ig^a'/3'(y + )^(y + ) + ^V2a'/3(y+)</>" (y+) 
+ |V<Ky + )B + -V) + |^«'5 + - / (y + )^(y + ) - ^|a'5+-(y+)^(y + ) 



+ ^B-~(y+)<p(y + ) + ^V2^ 3 S+ (3) (y+)0(y+) 

- |LvWV V¥(v + ) + ^vWiJ+V) W 

256 256 

- —V2^(y + )B-'(y + ) + —V2V^B-(y+)<f>'(y + ) 



16 



27 



^/2a>(y+)0"(y+) + J^aW) 2 } 
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+i^+{||a' 2 5+- (3 VM y +) - ^B+-'(yW) 
+ |^ + -(v + )* (3) (v + ) + ^4>{y+)B-\y+) 
~ ^B-(y^'{y^) + ^V2^ 3 Hy + )B-"(y" 



^V^V V^V) + ^V^V W(y + )} 

+^ 3 (^ + ) 2 {^« /2 0(y + )5- ,/ ( y + ) - ^|a' 2 s--'(^)^( y +) 

' 64 , 2 . 

r 

729 



176 
729 

64 _.___„. 12S 

1 (y ) - 7 

+ ^-«- J B-(y + )^(y+)} 
Equation of motion for (12 terms): 

a' 

QIC Kfi c _ 

- ^V^V^G/ + )£ + V) + ^V2V^£ + (y + )0V) + 
+^ + {§^>(y+)£+- V) - §W-(^'( y +) 
-^V^B-(y + )#y + )} 



Equation of motion for F (11 terms): 



+ — V2V^B-(y+)<P(y+)} 
+ K 3 (V + ) 2 {^a>B-(y+)4>(y + )} 



Equation of motion for (3 (11 terms): 

= - V + p'(x+) - 
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+ K*{^V2F(y + )^y + ) + ^/3(y+)<Xy+) - —V2B+-(y+)^y 



+ 



^^4>{y + )B+'{y+)- l ±<P{y+f} 



I 729 



o 



-)B+- / (y+) + ^V2a>B+-(y+)4>'(y' 



+ ^^B-(y + )4>{y+)} 
+ KHV+f{^V2 a 'B-(y+My+)} 



Equation of motion for B + (20 terms): 

B~(x + ) 



0= -V + B'\x + ) 



a 1 



+K 3 { - -^-V2^(y + )F>(y+) - -^-V2V^F(y+)0>(y+) - J±v/^>(y+)/3' \y + ) 



_V^/3(y+)^(y+) + ^V2V^(y+)£+- V) 
572 



+ —V2^B + -{y+W(y + ) - JL a >ct>(y + )B + " (y+) - -a* B + ' \y+)<f>' \y+) 
256 



- ^B-(y+)<P{y + ) + g^(y + y(y+)} 
64 



32 

81 ' 
256 



943 -V2V^ + — yftV* B" {y + )<j>{y + ) 

44/44 i 
- — aV>(y + )£ - —a'B-{y+)<t>\y+)) 

+ K 3 (V+) 2 {-^V2^ 3 Hy + )B-'(y + ) - ^V2^ 3 B-(y+W(y+)} 



Equation of motion for B (12 terms): 

B+{x+) 







V + B + (x- 



{-w?, B+{v+)<t,{v+) } 



20 



+K ? >V + {^V2^F(y+)4>{y + ) + ' f3(y + )cf>(y + ) + '^V2^/B+- (y+)^ 

5 



+ ^'4>{y + )B + \y + ) - ^V^cf>(y + ) 2 } 



+K\V+f{^V2^ 3 Hy + )B + -'(y + ) + i|V2^ 3 i? + -(y + )^( y + ) 



110 

H a 

729 



/ 5-(y + )0(y+)} 



+ ^ 3 (^ + ) 3 {^V2^ 3 S-(y+)0(y+)} 
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D. Determinant of polynomial matrices 



In this paper, we have to evaluate the determinant of a matrix M whose entries are poly- 
nomials in two variables u and y. The result is obviously a polynomial in u and y as well. 
One could of course use the traditional methods, for example LU decomposition, but one 
then notices that the complexity of the calculation grows very fast with the matrix size. 
This is easy to understand when one notes that such methods require the matrix to be over 
a field. One will then have to think of elements of M as rational functions of u) and y. The 
calculation of the determinant will then involve potentially heavy computations on large 
fractions. Moreover, since the determinant is a polynomial, the end result will involve a 
large simplification between numerator and denominator, meaning that we might be doing 
many more calculations than necessary. In fact we observed that, while mathematica can 
calculate such a determinant for a seven by seven matrix, it fails to do so for a fifty by fifty 
matrix whose entries are polynomials of degrees at most eight in u> and at most one in y. 
This prompted us to look for alternative methods. 



An elegant method [48| is based on the discrete Fourier transform. We mention also 
that there is another method |19| , that can efficiently calculate the degree of the determi- 
nant. This is in fact what we need, but it is unclear to us how to generalize this method 
to polynomial matrices in several variables, so we will use the first method based on the 
Fourier transform. 

To explain this method, it will be enough to concentrate on the case of polynomials 
in one variable x; the generalization to several variables will then be obvious. The crucial 
fact to realize is that polynomial multiplication is, in some sense, a convolution. More 
explicitly, let's consider two polynomials p{x) and q(x) of degrees d± and e?2 

rfi d% 
p(x) = J2PnX n and q(x) = ^ q m x m . (D.l) 

n=0 m=0 

We define N = d% + d% + 1 and we write the coefficients p n and q n into vectors p and q, 
both of length N. 

P = {po,Pi,---,PN-i) , q= (qo,qi,...,qN-i), (D.2) 

where p n = for n > d\ and q m = for m > d<i- We can then write the product t(x) = 
p(x)q(x) also in a vector of the same length because it has degree N — 1. The components 
of t are then given by t n = Ylm=o PmQn-m- If we define the negative components of p and 
q by imposing the periodicity p n = pn+u (and similarly for q n and t n ), we can extend the 
summation to the whole range 

N-l 

t n='^2 PrnQn-m- (D.3) 
m=0 

This in fact doesn't change the sum because all additional terms are zero; but in the form 
(|D.3| ), one recognizes that t is the cyclic convolution of p and q 

t = p * q. (D.4) 



- 41 - 



We can now use the fact that the discrete Fourier transform defined by 

JV-l 

T( a ) = (A ,A 1 ,...,A N - 1 ), with A k = £ e 27Tikn/N a n . (D.5) 

n=0 

changes cyclic convolution into multiplication. Namely 

F(t) = J-(p)J-(q), (D.6) 

where the multiplication of two vectors is done component-wise. One can then perform an 
inverse discrete Fourier transform to obtain t(x). 

Let us now consider, for a little while, the Leibniz formula for the determinant d(x) of 
the square matrix M(x) of size n 

n 

d(x) = sgn(^) II M ^d) i x ) • (D-7) 

We estimate an upper bound N — 1 on the degree of d(x) and write its coefficients in the 
vector d of length N; and we also write the coefficients of Mi Ax) in vectors Mj each of 
length N. We can now write the polynomial product in (p.7|) as a usual product of the 
components of the discrete Fourier transforms of Mjj 

n 

F{d) = sgn(a)l\T(M iMi) ). (D.8) 

o-eSVi i=l 

Now if we define the matrices M n by 

(M n ) hJ = (-F(M tJ ))n, (D.9) 

the Leibniz formula becomes 

n 

(.F(d)) n = ^ sgn( C r) n^nk^)- (D.10) 

In other words, the components of ^(d) are the determinants of the matrices A4 n whose en- 
tries are numbers; these determinants can therefore be calculated with a standard method. 
To summarize, the algorithm to calculate the determinant d(x) of M{x) is: 

1. Find an upper bound iV — 1 on the degree of d(x). 

2. Calculate the matrices A4 n for n = 0, . . . , N — 1, whose entries are the components 
n of the discrete Fourier transform of the entries of M. 

3. Calculate the determinants of the matrices A4 n . These are the components of J-(d). 

4. Inverse discrete Fourier transform J-(d) in order to find the coefficients of the poly- 
nomial d(x). 
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We make a few comments. First, it is well known that the discrete Fourier transform 
can be made extremely fast with the Fast Fourier Transform algorithm, so this method 
is relatively efficient. On the other hand, if one is interested in determining the degree of 
d(x), numerical errors can make it hard to decide whether a term ex n appearing in the 
answer for the determinant, with e a very small number in absolute value, should be kept 
or whether it comes from numerical imprecision. For this reason, it might be necessary 
to run the algorithm with a precision of many digits. At last, this algorithm generalizes 
almost trivially to polynomials in several variables; the discrete Fourier transforms have 
then to be replaced by multidimensional discrete Fourier transforms. 

E. Review of linear dilaton CFT 

In the following we want to review briefly how the linear dilaton CFT arises and quote the 
most useful results that are needed to compute correlators and the equations of motion in 
a linear dilaton background. A similar, even shorter overview can be found in Section 
2. Much more is covered in [|50[1 , Sections 2.5, 2.7 and 3.7, where additionally the role of 
the dilaton field in string interactions is illuminated. 

Let's start by considering the description of string theory in curved spacetime. After 
all string theory is supposed to give a framework for quantum gravity. As a reminder the 
Polyakov action that serves as the starting point for string theory is 

Sp = -^-f I d 2 aV=7 7 ab d a X»d b X»^ u . 
4vra' J M 

As a natural extension one is tempted to replace the flat Minkowski metric r]^ u with a 
general metric G^ u (X). Considering spacetime as a coherent background of gravitons is 
slightly dubious, as it means string theory is formulated in a background of strings. But 
let's try anyway. If we add the graviton background we might as well add the background 
from the two other closed string massless states, the antisymmetric B^ v tensor and the 
dilaton resulting in a so called nonlinear sigma model action 

Sa = / d2<7 ^ { (j ab °^ w + ieabJ3 ^ (*)) + a ' m w} . 

where 7 ab , R are the worldsheet metric and curvature. 

If we now demand this action to be Weyl invariant also at the quantum level the 
corresponding /3-functions for G,B,& of the renormalization group flow have to vanish 
(<^scale invariance). To first order in a', only derivatives of the dilaton $ appear in the 
/3-functions. This means changing the dilaton by a constant is allowed. A constant shift 
in $ yields a contribution proportional to the worldsheet Euler number x> defined as 

X = --r^- 1 [ d 2 aV=jR- 
4vra' J M 

In fact the Euler number determines the interaction strength, such that the open string 
coupling constant g is related to the expectation value of the dilaton by 
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Typically in quantum field theory the coupling is a free parameter, so different values for the 
coupling represent different theories in sharp contrast to string theory: different coupling is 
just a different background of the same theory. However since we cannot determine these 
background values from the dynamics, this distinction moves the difficulty elsewhere, but 
in principle, there is no free parameter as in QFT. 

If the dilaton is not constant, then in order to preserve Weyl invariance, the critical 
dimension of the string theory may be altered. One choice of fields that preserves Weyl 
invariance is 

G llv {X) = Vlxv , B^ = 0, ^{X) = V ll X^ 

where $ is linear in X (linear dilaton), and the constant vector is the dilaton gradient. 
It singles out a particular direction in spacetime, breaking spacetime translation symmetry. 
Momentum conservation is then modified as 

3 3 

The (3 functions vanish if 

Conversely the central charge of this theory is 

c = D + 6d%V. 

Hence different values of the dilaton gradient V give different CFTs. We will choose a 
light-like dilaton gradient, VpV^ = 0, then the central charge and critical dimension are 
unaffected. 

The complete (Euclidean) action for the linear dilaton on worldsheet S with boundary 
<9X then becomes 

SwS = i^ 7 L d2fJ ^ ( 7 abl l^aX»d b X u + ctRVpX*) +^J ds dskV^, (E.2) 

where k denotes the worldsheet boundary curvature. Varying the action with respect to 
the metric gives the energy-momentum tensor 

T (z) = -1 : dX^dX* : +V^X", (E.3) 

or 

the usual expression with an additional dependent contribution. 

Gauge fixing the metric permits obtaining the same metric as before on the complex 
plane, ds 2 = dzdz such that the gauge fixed action has the same form as in the free theory 
with = 0: 

This theory deserves the name CFT, for the action ([E.4D is invariant under the conformal 
change of coordinate z — > f(z), where the field transforms as 



X^ 



(z, ^/ol" (z, z) = X» (/ (z) , / (*)) + |V log \f'(z)\ 2 . (E.5) 
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Note that this transformation has again an extra piece containing the dilaton gradient. We 
want the (real) boundary of the upper half plane to be mapped into itself so we can use 
the doubling trick and (z,z) — > X^ (z), thus we require f (z) = f(z). Denoting the 
real coordinate with y G M in contrast to complex z, a boundary field transforms as 

X» (y)^X^ (f(y)) + a'Vnog\f (y)\ , 

with the visible difference that the 2s have disappeared. Expanding the energy-momentum 
tensor into Laurent modes, we find the new (matter) Virasoros 



1 °° / ^ 



(E.6) 



The conformal weights of any operator can be determined from the OPE with the energy- 
momentum tensor. For example, the weight of : e tk ' x ^ z ' z ^ : is 

,(k 2 

a — + 1— ^ — 



On the boundary, things are a little more complicated. The normal ordering 

: X» {z u zt) X v (z 2 , z 2 ) : = X» {z u z{) X v (z 2 , z 2 ) + |V log \zi - z 2 \ 2 

is just valid in the bulk, not on the boundary. On the real line the Green's function 
Gi2 = -%- log I z\ — z 2 \ needs to be modified (call it G' 12 ) to satisfy the Neumann boundary 
conditions (<9 — d) G' l2 = 0. This is accomplished by an additional image charge term 

G'xi = -y log \zi - z 2 \ 2 - y log \zi - z 2 \ 2 . (E.7) 

Normal ordering should remove the divergence of this Green's function at zero separation, 
but the usual subtraction ^-r/^" log \z\ — z 2 \ is not enough, it only cancels the first diver- 
gent term, thus it has to be doubled. The new prescription goes by the name of boundary 
normal ordering, distinguished from the bulk conformal normal ordering by new symbols 

IX» ( Vl ) X v (y 2 ) I = X» ( Vl ) X v (y 2 ) + 2a' log \ Vl - y 2 \ . 

We stress again that the reason for introducing another normal ordering is just to have 
well-behaved (finite) expectation values for products of boundary normal ordered operators. 

When taking real derivatives of boundary operators, factors of 2 need to be taken into 
account because d z = \d y . This is a source of great confusion. Furthermore operators 
on the boundary may also have different weights from their "relatives" defined on the 
interior/bulk because of that, e.g. * e lk ' x ( z )* has weight 



a 



'(k 2 + ik-V). (E.8) 



It is instructive to check the consistency: the linear dilaton is a background field, thus 
non-dynamical. If we set = in the above formulae the free string theory is recovered, 
as it should. 
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